Renormalization of the upper critical field by superconducting fluctuations 
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We study the effect of superconducting fluctuations on the upper critical field of a disordered 
superconducting film at low temperatures. The first order fluctuation correction is found explicitly. 
In the framework of the perturbative analysis, superconducting fluctuations are shown to shift 
the upper critical field line toward lower fields and do not lead to an upward curvature. Higher 
order corrections to the quadratic term coefficient in the Ginzburg-Landau free energy functional 
are studied. We extract a family of the mostly divergent diagrams and formulate a general rule 
of calculating a diagram of an arbitrary order. We find that the singularity gets more severe with 
increasing perturbation theory order. We conclude that the renormalization of the Ginzburg-Landau 
coefficients by superconducting fluctuations is an essentially non-perturbative effect. As a result, 
the genuine transition line may be strongly shifted from the classical mean-field curve in a two- 
dimensional superconductor. 
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I. INTRODUCTION 

The magnetic field-temperature phase diagram of type 
II superconductors has been of experimental and theo- 
retical interest for a very long time. Much interest has 
recently focused on the disordered two-dimensional su- 
perconductors partly motivated by the two-dimensional 
nature of the high-T c cuprates although disordered thin 
films of "ordinary" superconductors have been studied for 
years. The subject matter of this paper is the tempera- 
ture dependence of the upper critical field, H C 2(T), above 
which the system goes normal. In particular, we consider 
the problem of superconducting fluctuation effects on the 
upper critical field of a disordered superconducting film 
at low temperatures. 

There have been a number of experiments studying the 
behavior of the upper critical field line as a function of 
temperature and disorder. Some of themi have shown 
strong deviations from the classical mean-field theory^ 
such as an upward curvature in the H C 2(T)-\me. These 
experiments have induced a considerable theoretical in- 
terest in the subject. 

Let us start with recalling the previous theoretical 
studies. The first attempt to explain the observed effects 
has been made by Golubov and Dorin£ who considered 
the combined effects of disorder and Coulomb interaction 
in first-order perturbation theory. Golubov and Dorin's 
calculations predicted some upward curvature in the up- 
per critical field line. However, later re-examination of 
the Golubov and Dorin's theory showed some very tech- 
nical deficiencies in the original calculations. Namely, 
Smith et alA have performed both first-order perturba- 
tion theory calculations and Oreg and Finkclshtcin renor- 
malization group treatment of the problem. Their care- 
ful analysis has not shown any anomalous upward curva- 
ture in either case. So, the final verdict was that disorder 
and Coulomb interaction together could not explain the 
observed deviations from the mean field results. 

Zhou and Spivak have suggested 6 " a mechanism of H C 2 



enhancement based on a possibility of the spontaneous 
formation of a granular-like superconducting structure in 
a disordered film at a low enough temperature. The idea 
is that optimal fluctuations in the distribution of impuri- 
ties may lead to the formation of local regions where the 
local critical field exceeds the system-wide average value. 
At a low enough temperature, these disorder-induced su- 
perconducting droplets form a Josephson network allow- 
ing for the global superconductivity in the film. Recently, 
Galitski and Larking have included the effects of quan- 
tum fluctuations into the picture and showed that the 
"mesoscopic disorder" (i.e., usual non-magnetic impu- 
rities) was very weak to produce a considerable effect 
on the critical field line and the global superconducting 
glassy state suggested by Zhou and Spivak 6 should be de- 
stroyed by quantum fluctuations. However, the "pinning 
disorder" (grain boundaries, dislocation clusters, etc.) 
may be strong enough to compete with the quantum 
fluctuations. The prediction is that the anomalous up- 
ward curvature due to the optimal disorder fluctuations 
is possible but the corresponding effect is non-universal 
and depends on the pinning properties of the system. 
Let us mention a Commem™ on the paper— by Ikeda, 
who emphasized a possible importance of superconduct- 
ing fluctuations. The purpose of the present work, which 
is partially motivated by Rcf. 8], is to explicitly evaluate 
fluctuation corrections to the upper critical field line at 
low temperatures. 

Let us emphasize that the definition of the supercon- 
ducting transition point itself is a very delicate issue in 
a two-dimensional system. Strictly speaking, at finite 
fields, superconductivity as a state with zero resistivity 
is never achieved which is basically a manifestation of the 
general Mcrmin- Wagner theorem^ In real experiments, 
what is observed is a drop of the resistance to some very 
small but finite values. Such a point is usually defined 
to be the transition point. Further decrease of the tem- 
perature or external field yields a very slow decay of the 
resisitivity. From the theoretical point of view, the def- 
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inition of the critical point is also non-trivial. One can 
define the transition point as a point at which the super- 
fluid density p s becomes non-zero. Another approach is 
to study the behavior of the quantity (ip(O)ip(r)} = A(r) 
by the use of Ginzburg-Landau type expansions and de- 
fine the transition point as a point at which the quadratic 
term vanishes. There are no solid grounds to believe that 
these two approaches are equivalent. In the present pa- 
per, we follow the latter approach, which, as we think, 
should correspond to the point at which the critical drop 
in the resistance to some finite value takes place. 

The question of T c renormalization in zero field has 
been recently addressed by Larkin and Varlamov in the 
Review The two-dimensional result for the shift 

of the transition temperature was found as 5T C /T C ~ 
[Gi InGi], where Gi is the Ginzburg parameter, which in 
a disordered two-dimensional superconductor is of order 
of the inverse conductance. Thus, one concludes that the 
effect of fluctuations in zero field is quite small and the 
fluctuations can move the transition point only within an 
extremely narrow Ginzburg region. We will return to the 
discussion of the validity of this conclusion at the end of 
the paper. 

Our paper is structured as follows: In Sec. II, we for- 
mulate the Ginzburg-Landau theory for the case of a 
strong external magnetic field. We emphasize that the 
Ginzburg-Landau expansion in this case can be made 
with respect to the modulus of the order parameter only 
and not on its spatial gradients. However, one can still 
write down a formal operator expansion. The vanishing 
of the lowest Landau level matrix element of the oper- 
ator A appearing in the quadratic term, corresponds to 
the transition point. In Sec. Ill, we study the effect of 
superconducting fluctuations on the upper critical field 
in first order perturbation theory. We calculate the first 
correction explicitly. We find that the first correction 
alone can not lead to a monotonous upward curvature of 
the H C 2(T) line but can yield a decrease of the critical 
field compared to the mean-field valued In Sec. IV, we 
study the higher order perturbation theory corrections. 
We extract a family of the mostly divergent diagrams 
the effect of which parametrically exceeds the first and 
second order contributions. The first term in the fam- 
ily appears in third order perturbation theory only. We 
show that the singularity gets progressively stronger in 
higher orders. We conclude that the renormalization of 
the Ginzburg-Landau coefficients and, consequently, the 
transition line by superconducting fluctuations is an es- 
sentially non-perturbative effect. 



II. GINZBURG-LANDAU EXPANSION IN 
STRONG FIELDS 

It is well-known that a superconductor in a relatively 
weak magnetic field can be described with the aid of the 
Ginzburg-Landau theory. The Ginzburg-Landau expan- 



sion of the free energy has the following standard form: 
T[A(r)]=T N + |d 2 r|a|A(r)| 2 + 6|A(r)| 4 (1) 

+ ^-|[V-2z e A(r)]A(r)| 2 j, 

where J-n is the free energy of a normal metal, A(r) is 
the superconducting order parameter, a, b, and m are 
some coefficients which depend both on temperature and 
magnetic field. Coefficient a vanishes at the transition 
point. 

To derive Eq.([Q) from the microscopic theory, i.e. from 
the initial BCS hamiltonian, one has to fulfill the follow- 
ing standard steps: First, the quartic interaction term 
should be decoupled with a Hubbard-Stratonovich field 
A. Then, the functional integral on the electron degrees 
of freedom becomes Gaussian and can be easily evalu- 
ated. The resulting action 5(A) is quite complicated, 
however in the vicinity of the superconducting transition, 
A is small and the action can be expanded on it. If one 
is interested in the equilibrium thermodynamic proper- 
ties of a superconductor, the time-dependence of A may 
be suppressed and one can talk of deriving an effective 
free energy functional rather than of an effective action. 
However, another essential assumption has to be made in 
order to get Eq.lJIJ. Namely, one has to suppose that the 
spatial variations of the order parameter are small enough 
so that term |(V — 2ieA(r)) A(r)| is small in some sense. 
Finally, expanding both on the order parameter and on 
its spatial gradients on can derive the Ginzburg-Landau 
free energy functional Q the minimum of which yields 
the famous Ginzburg-Landau equations. 

The above arguing partially fails in strong fields 
when the spatial dependence of the order parameter be- 
comes important. In strong fields, the expansion on 
|(V — 2ieA(r)) A(r)| is no longer possible. However, one 
can still expand on the magnitude of the order param- 
eter. The corresponding saddle point equation can be 
written as: 

iA(r)=|c(r,r')A(r')dV (2) 

- j B(r, r 2 , r 3 , r 4 )A*(r 2 )A(r3)A*(r 4 )d 2 r 2 d 2 r 3 d 2 r4. 

Here C is Cooperon (see Fig. I) which is a linear operator 
in the case under consideration. In the presence of a 
magnetic field the operator can be written as 

C(r,r')=Tj2Qe(r)G- s (r'), 

e 

where Q is the electron Green function in the magnetic 
field and e is the fermion Matsubara frequency. 

Non-linear operator B is described by the square dia- 
grams similar to the ones shown in Fig. 2. Let us note, 
that this operator is not singular at the transition point 
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FIG. 1: Diagrammatic equation for the fluctuation propagator © (curly line). 
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FIG. 2: First order fluctuation corrections to the quadratic part of the Ginzburg-Landau free energy expansion. 



and in its vicinity and can be considered as a local quan- 
tity for a dirty superconductor [see explicit expression 
iJTS) below]. 

In the semiclassical approximation, the magnetic field 
dependence in the Cooperon can be factored out: 



C(r, r') = C* (0) (r - r') cxp I -2ie J A(s)ds 



(3) 



where (r) is the Cooperon without magnetic field. 
Formally, operator C can be written in the following 
form: 



C = J C (0) (r) exp{-ir7r}d 2 r, 



(4) 



where 7r = — zV — 2eA(f) is the operator of the kinetic 
momentum, which can be expressed in terms of the cre- 
ation and annihilation operators in the Landau basis. 
Again, in small fields, one can expand the exponent in 
Eq. J3J which will result in the gradient term ir 2 [see 
Eq.©]. In strong fields, we can not do such an expan- 
sion. 

The transition point is defined as a point at which 
the coefficient in front of the A 2 term vanishes. The 
corresponding operator has the form: 



A = X~ 1 -C. 



(5) 



The inverse operator can be easily recognized as the pair- 
ing vertex which corresponds to the ladder summation as 
shown in Fig. 1. 



£(0) = k 



(6) 



The divergence of this quantity at the transition point 
corresponds to the BCS instability. 

It is possible to calculate pairing vertex £(f2), which is 
also called the fluctuation propagator, as a function of the 
total energy in the Cooper channel £1 at arbitrary mag- 
netic fields. From Eq. JI}, one can see that the Cooperon 



can be written as a series containing even powers of ki- 
netic momentum 7r only. Thus, the Cooperon is a di- 
agonal operator in the Landau basis, which makes [see 
Eqs. il -3 1( > fl J the fluctuation propagator diagonal as well. 
The corresponding matrix elements have the following 
explicit form: 



1 



JV(0) 



hi- 



T 



c0 
f 

2 



101 



4ttT 



(7) 



where N(0) is the density of states at the Fermi-line, T c o 
is the transition temperature in zero field, Q = 2imT 
is the Matsubara frequency, and = 4eT>H, with T> 
being the diffusion coefficient and H the external mag- 
netic field. The upper critical field line is defined by the 
condition: 



A = £ - 1 (0) = 0, 



(8) 



where index "0" corresponds the lowest Landau level ma- 
trix element. The above equation yields the well-known 
Gor'kov's mean-field curve for H C 2(T) (see solid line in 
Fig. 3). 

In order to study the effect of superconducting fluctu- 
ations on the H C 2 transition line, one has to go beyond 
the mean field theory. The first fluctuation correction is 
described by the diagrams shown in Fig. 2. The corre- 
sponding correction to the Ginzburg-Landau coefficient 
is 

W (1) (n,rj) -T 2 ^B £ ,o(ri,r 2 ,r3,r 4 )/:n(r3,r4), (9) 

where B is the impurity box (see Fig. 2) and Cq is the 
fluctuation propagator introduced in the previous sec- 
tion. 

Let us note that operator 6 A in Eq. © is diagonal in 
the Landau representation, since it can be written in the 
form analogous to Eq.©. The fluctuation propagator 
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Temperature, T/T c , 



FIG. 3: Renormalization of the upper critical field by superconducting fluctuations. The solid line shows the mean-field 
Gor'kov-Helfand-Werthamer result 2 . The dashed line is the upper critical field with the account for the superconducting 
fluctuations in first order perturbation theory. 



can be expanded on the Landau wave- functions as follows 
(we use the Landau gauge here): 

f d 

£(ri,r 2 )=^£„ / ^VvOOVwfa)) (10) 



where 



Vv(r) = e l ™ ^ sc \x - x ), x =pL 2 H , (11) 



with ^n° sc '' (x) being the harmonic oscillator eigenfucntion 



osc )(<r\ — 



b£ sc >(x) 



v / 2 n nl 



1/4 



•^H„( — ), (12) 



and Lh = (2eH)~ 1 the magnetic length for the fluctuat- 
ing Cooper pairs. 

Thus, the first fluctuation correction can be written as: 



e.Q n 



r 

SA^ =T 2 J2J2 / II ^ ^,o(ri,r 2 ,r 3 ,r 4 ) 

dkdk' 

— — T tpnk(ri)^ nk (r 2 ) 
(2tt) 

x Vofe' (r 3 )ipo p (r i ) . (13) 

At the superconducting transition point, only the low- 
est Landau level matrix element of the fluctuation prop- 
agator is divergent. In the vicinity of the transition it 
has the simple form: 



1 \[H-H c2 (T)] 2|fi| 



AT(O) 



H c2 (0) 



H 



(14) 



where H c2 (T) is the mean- field value of the upper crit- 
ical field. Since only the lowest Landau level provides 
a divergent contribution, one can keep just one term in 
the sum over n in Eq. I|13|) . With the same accuracy, one 
can put 17 = in the expression for the box diagrams 
as the divergence comes from the small values of Q in 
the sum. In this approximation, the ladder summation 
in diagram "c" in Fig. 1 is not necessary and we are left 
with the Hikami-box like diagram, the result for which is 
known 



S £ n=o(ri,r 2 ;r 3 ,r 4 ) 



tt7V(0) 



n 



i 



x5(ri - r 2 )<5(ri - r 3 )<S(ri - r 4 ) 



where 



| e | + -2^[0 (1) -0 (3) ] + [0 (2) - 3(4)]' 



fl(fc) = -iV-2e(-l) fe A(r). 



(fe) 



, (15) 



\vd 2 



come from the 



In Ea. (|15|l . the factors 

four impurity vertices, three delta-functions correspond 
to the locality of the quantity in the dirty limit (tT c q <C 
1, where r is the scattering time), and the last factor 
comes from calculating the square diagrams themselves. 
Let us note that the leading frequency and momentum 
independent terms are cancelled out and one has to keep 
e and finite to get a non-zero contribution. This prop- 
erty can be proven not only for box-like diagrams shown 
in Fig. 2, but for any 2n-sided polygon. 

Acting by operator (|15fl on the lowest Landau level 
eigenfunctions in Eq.@ and calculating the Matsubara 
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FIG. 4: Second order fluctuation corrections to the quadratic part of the Ginzburg-Landau free energy expansion. 




FIG. 5: Some higher order diagrams are shown. Graphs "b" and "c" are examples of reducible diagrams which can be separated 
into several irreducible parts by cutting curly lines only. 



sum, we derive the following expression: 



SA 



(i) 



1 



2nL 2 H (4irT) 2 \8nT 



(16) 



The last factor can be evaluated with the logarithmic 
accuracy as follows: 



CI 



H 



4tt7V(0) 



In- 



n 



H 



H - H c2 (T) n H 



T T V H c2 (0) 

T H c2 {0) 
2fl H H- H c2 



T 



(17) 



Let us remember that N(0) is the density of states at 
the Fermi- line, Q H = AeVH = ^-T c0 , where V is the 
diffusion coefficient, e is the electron charge, H is the 
external field, T c q is the transition temperature in zero 
field, and 7 ss 1.78 is Euler's constant. H c2 (T) stands 
for the mean field temperature dependence of the upper 
critical field which is defined by Eas. i|7l8|) and in the low- 
temperature limit has the form: 



H c2 {T) « H c2 (0) 



1 - 



27 2 



T 



Eas. (|16ll7|l determine the first order correction to the 
upper critical field which is defined by the following equa- 
tion: 



^0 + SA y 



(i) 



0. 



This is a non-linear algebraic equation, which can be eas- 
ily solved by numerical means. As one can see in Fig. 3, 
the first correction yields a lower upper critical field com- 
pared to the mean field result^ As expected^ the per- 
turbative analysis does not lead to an upward curvature 
in the ff c2 (T)-line. 

Let us study the asymptotic behavior of the ob- 
tained expressions. If the following inequality is satisfied 
T/T c0 < [H - H c2 (T)} /H c2 (0) < 1, one can obtain the 
following zero-temperature result: 



^ = Gi ln 
N(0) 4tt 



H 



H-H c2 (0) 



(18) 



where we have introduced the dirty limit Ginzburg pa- 
rameter as Gi = (eft) -1 -c 1. One can see that the 
first order correction is only logarithmic. In the case of 
small but finite temperatures [H — H c2 (T)] / H c2 (0) 
T/T c0 < 1 we have 



SA 



(i) 



N(0) 



7G1 
T2" 



H 



H-H c2 {T) 



T 

T c q 



(19) 



We see that in this temperature domain, the effect of the 
superconducting fluctuations is more pronounced. 



III. HIGHER ORDER CORRECTIONS 

In this section, we will study higher order corrections 
to the mean field result. For the sake of technical sim- 
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plicity, we will focus on the zero-temperature case which 
we will use to estimate the divergence rate of higher or- 
der diagrams. Generalization to non-zero temperatures 
is straightforward. 

The simplest way to construct a higher order correc- 
tion is to study a diagram made of a polygon with all 
possible impurity avergaings and n curly lines connect- 
ing 2n vertices of the (2n + 2)-sided polygon. In Fig. 4, 
all possible topologically non-equivalent second order di- 
agrams are shown (see also Fig. 5, diagram "a"). Note 
that in the highly disordered limit, these polygons are 
local quantities which can be written in the form sim- 
ilar to (fH)l) . The contributions due to such diagrams 
can be easily estimated. Expressions for these diagrams 
contain n integrals over the frequencies running through 
n fluctuation vertices, which leads to a contribution of 



order of In™ 



-Hc2 



Therefore, we conclude that per- 



h-h c2 

turbation series due to the single-polygon (irreducible) 



where n is 



graphs has the form J2 n C"GIi n ln r 

the perturbation series order and C n are some combi- 
natorial factors. Let us mention that similar diagrams 
have been studied by Kee et alMi who constructed a 
non-perturbative resummation technique to sum up the 
dominant contributions. In our problem, such a resum- 
mation is not possible since, as we shall see, another type 
of diagrams delivers the dominant contributions. 

Problems start to appear in third order perturbation 
theory when one is able to construct a diagram contain- 
ing several polygons connected by curly lines only (see, 
e.g., Fig. 4; diagram "b"). We will call a diagram re- 
ducible if it can be separated into several parts by cut- 
ting curly lines only. Otherwise, we will call a diagram 
irreducible. Let us consider a general case of a graph in 
n order perturbation theory which contains r irreducible 
parts. In this case, simple dimensional analysis of the 
problem yields the conclusion that we have (n + 1 — 7') 
integrals running through n singular fluctuation propaga- 
tors. Therefore, the contribution due to such a diagram 
can be estimated as follows: 



5A, 



(n,r) 



N(0) 



Gi r ' 



H-H c2 



(r-l) 



In 



(n- 



H-H, 



One can easily see that in higher orders the singular- 
ity gets enhanced and contains a power low dependence 
on the closeness to the transition compared to only 
logarithmic dependence l|18|) in first and second order 
perturbation theory. At low but finite temperatures 
[H - H c2 (T)} /H c2 (0) < T/T c0 < 1, we obtain the fol- 
lowing estimate for the contribution of an n order graph 
built up of r irreducible fragments: 



SA 



(n,r) 



N(0) 



(n-r+1) 



H c2 (Q) 



H-H c2 {T) 



It is clear that the more irreducible parts contains a di- 
agram, the more singular contribution it derives. There- 
fore, mostly divergent diagrams should contain as many 



cubic irreducible elements as possible. The correspond- 
ing subseries is built up of chain and ring-like diagrams 
(see Fig. 4, diagrams "b" and "c"). The contribution 
from a diagram which contains r Hikami boxes is found 
as (T = 0): 



SA 



(2r-l,r) 

N(0) 



H c2 


r-l 


H — H c2 





H — H c2 

Such diagrams provide dominant contributions in the 
expansion. From Ea.l j20|l . we see that although the 
singularity gets progressively stronger in the higher or- 
ders, it tends to a finite limit. Namely one can 
conclude (by taking the limit n — > oo) that if 
[(H-H c2 )/H c2 ]ln- ¥ [H c2 /(H-H c2 )} » Gi 2 , each 
term in the perturbation series is small. 

Let us mention that keeping any finite number of terms 
in the perturbation series does not make any qualitative 
changes to the shape of the renormalized H c2 (T)-cmve 
compared to the first correction and to the mean-field 
results (see Fig. 3). However, increasing the number of 
terms in the series leads to a further decrease of the up- 
per critical field. This gives a very tentative indication 
that superconducting fluctuations in a two-dimensional 
system can significantly change the upper critical field 
but, probably, should not yield an upward curvature. 



IV. CONCLUSION 

We have considered the effect of superconducting fluc- 
tuations on the upper critical field of a two-dimensional 
disordered superconductor at low temperatures. We have 
performed an exact calculation of the first perturbation 
correction and found the corresponding shift of the crit- 
ical line. We have also estimated higher order correc- 
tions. The renormalized upper critical field was shown 
to be shifted toward lower fields compared to the classi- 
cal mean field curve. However, no upward curvature was 
found in the framework of any finite-order perturation 
theory expansion. 

We have carefully studied higher order contributions 
and found that each consecutive order in the perturba- 
tion theory enhances the singularity. We formulated a 
general rule of constructing and calculating the domi- 
nant contributions in the perturbation series. One of our 
main conclusions is that the derivation of the Ginzburg- 
Landau functional in the case of strong fields is hardly 
possible in the framework of the conventional perturba- 
tion theory technique. 

Indeed, our results and conclusions are directly appli- 
cable for the case of low temperatures and strong fields 
only. However, we feel that the issue of the perturbation 
theory applicability and convergence may be essential in 
the case of weak or even zero field as well. Certainly, in 
this case calculations are technically absolutely different 
from the case of strong fields, since one has to evaluate 
integrals over momenta instead of tracing over the Lan- 
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dau level indices. However, simple estimates show that 
reducible graphs bring up more singular contributions 
compared to the irreducible ones just like in the low- 
temperature case. Moreover, another problem of treating 
short- wavelength singularities comes out^ 

Taking into account these considerations we conclude 
that the problem of the renormalization of Ginzburg- 
Landau coefficients by superconducting fluctuations and, 
therefore, the renormalization of the transition point it- 
self, are essentially non-perturbative problems in a two- 
dimensional case. This makes the Ginzburg-Landau ap- 
proach ill-justified. By saying this, we certainly do not 
challenge the commonly accepted and well-tested phe- 
nomenological form of the Ginzburg-Landau functional. 



However, it is possible that the coefficients in the expan- 
sion are quite different from the conventional mean-field 
prediction due to the non-perturbative renormalization 
by superconducting fluctuations. 
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